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ABSTRACT 

The Effects of Shock Wave Precursors Ahead of 
Hypersonic Entry Vehicles. (December 1990) 

Scott Alan Stanley, B.S., Texas A&M University 
Chair of Advisory Committee: Dr. Leland Carlson 

A computational method and computer program have been developed to model 
the nonequilibrium shock wave precursor in a nitrogen gas. This method includes the 
effects of both chemical and thermal nonequilibrium. As part of this research, 
expressions have been derived for the mass production of each species due to 
photochemical reactions as well as for the effects of the absorption and emission of 
radiation on each individual energy mode of the gas. The shock layer and the radiative 
transfer solutions have been obtained by utilizing existing computational methods. The 
test cases presented represent entry conditions into the earth’s atmosphere for an 
aerocapture vehicle returning from Mars. 
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NOMENCLATURE 


AF - Radiation attenuation factor (-) 

£>; - Dissociation energy for the nth species (eV) 

e - Energy per unit mass (erg/g) 

e" - Electron/electronic energy (erg/g) 

E - Energy per particle (eV) 

E 2 - Second exponential integral (-) 
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h - Static enthalpy (erg/g) 
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k - Boltzmann’s Constant (1.38X10" 16 erg/°K) 
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n BF - Number of bound-free radiative processes 
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n BFi on - Number of bound-free photoionization processes 
- Number of molecular bands 
n s - Number of species 



Number density of the nth species (1/cm 3 ) 

Pressure (dyn/cm 2 ) 

Radiative flux (W/cm 2 sec) 

Universal gas constant (8.317xl0 7 erg/°K Mole) 

Heavy particle temperature (°K) 

Electron/electronic temperature (°K) 

Velocity in the precursor (cm/ sec) 

Mass production rate (g/cm 3 sec) 

Spatial variable in the precursor (cm) 

Absorption coefficient ratio as defined in Appendix B (-) 

One-half of the angle subtended by the body 
as defined in Figure 58 (degrees) 

Frequency (1/sec) 

Density (g/cm3) 

Optical depth (-) 

Change in the radiative flux affecting the 
energy mode indicated by M 

Change in the radiative flux due to absorption or 
emission through the radiative process indicated by P 



Subscripts 


e° - Zero point 

elct - Electronic 

i - For the ith process 

j - For the jth electronic level 

k - At the kth spatial point 

ke e . - Free electron kinetic 

n - For the nth species 

rot - Rotational 

tr - Translational 

vib - Vibrational 

oo - Freestream values 

v - At the frequency v 
Superscripts 

BB A4ina - Bound-bound atomic lines 
BB M harid - Bound-bound molecular band 
BF diss - Bound-free photodissociation process 
BF ion - Bound-free photoionization process 

FF - Free-free radiative process 

p - For a specific radiative process 
s - For a specific species 
TS - Using the tangent slab approximation 
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INTRODUCTION 

The recent emphasis placed on a manned mission to Mars and the return of men 
and samples to earth has caused an increased interest in the development of accurate 
methods for predicting the fluid flow around hypersonic entry vehicles. This interest is 
a result of the plan to use an aerocapture technique to provide the reduction in velocity 
necessary to place the spacecraft in earth orbit. The benefit of this technique is that 
aerodynamic drag, resulting from the interaction of the spacecraft with the earth’s 
atmosphere, instead of propulsive braking is used to slow the vehicle to orbital speeds. 

Such an approach provides a reduction in the fuel necessary for the mission and an 
increase in the return payload capabilities. The velocities experienced by a vehicle on 
return to earth from Mars are in the high hypersonic range, on the order of 1 1 Km/sec 
to 16 Km/sec. 1,2 

In recent years, the majority of work associated with Mars return has involved 
the shock layer, the region of hot gas behind the shock wave; an area which has received 
little attention is the shock wave precursor. The precursor is the region ahead of the 
shock wave in which radiation, primarily ultraviolet, emitted by the hot shock layer is 
reabsorbed by the gas. This absorption of radiation causes a heating of the gas in the 
precursor and the production of atoms as well as ions through photoionization and 
photodissociation reactions; these changes might also in turn affect the gas behind the 
shock wave. For example the preheating of the gas in the precursor as well as the 
introduction of electrons and ions could potentially increase the rate at which the gas 
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behind the shock approaches equilibrium; it has also been shown that for certain 
conditions the absorption of radiation ahead of the shock can cause significant increases 
in the radiative heating to the body. 3,4 Further, the presence of free electrons in the 
precursor also significantly affects communication with and identification of entry 
vehicles. 5,6 

The majority of the early work on precursors was performed using shock tubes 
and shock tunnels 5,7,8,9,10 ; in this early work, there was considerable debate as to 
whether the electrons observed ahead of the shock wave were due to diffusion through 
the shock front, photoionization, or tube wall effects. However, later research by 
Lederman and Wilson 6 as well as by Zivanovic 11 showed that the dominant cause of 
these electrons was indeed photoionization. In their work, Lederman and Wilson 
performed shock tube tests with a Venetian blind type insert in the tunnel which prevented 
the propagation of radiation but allowed the diffusion of electrons. With the insert in 
place, there were no electrons present ahead of the shock indicating that diffusion could 
not account for the presence of electrons in the precursor. Zivanovic performed tests in 
a ballistic range at velocities of up to 6.7 Km/sec, utilizing charged screens to prevent 
the diffusion of electrons; he found that these screens had no effect on the presence of 
electrons ahead of the shock. 

Another work of particular interest to this study was performed by Omura and 
Presley 12,13 . They conducted a shock tube study of the electron densities ahead of 
strong shock waves in nitrogen as well as in air mixtures. These studies by Omura and 
Presley are, to the authors knowledge, the only actual measurements made of the 
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precursor in air. In this work, Omura and Presley placed plastic filters in the shock tube 
to prevent the propagation of ultraviolet radiation. With these filters in place, there were 
no electrons present ahead of the shock wave; therefore showing that the precursor was 
due to photoionization of the gas resulting from the absorption of ultraviolet radiation. 

In the late sixties and early seventies, a number of analytical and computational 
studies were made in an attempt to model the precursor. Murty 14 performed a study 
that included the effects of the absorption of line radiation in the precursor in a hydrogen 
gas and Nelson 15,16 calculated the degree of ionization and excitation in the precursor 
ahead of a shock wave in an argon gas. Also, a number of studies including various 
effects were made by Clarke et al. for argon and hydrogen gases . 17,18,19 

Finally, the most recent work in shock wave precursors was by Tiwari and 
Szema . 20,21 In this work, they calculated the effect of the precursor on the shock layer 
and the radiative heating on a body entering into the atmosphere of Jupiter, which is 
composed of hydrogen. 

The purpose of the research described herein was to develop a technique for 
predicting the character and magnitude of the shock wave precursor ahead of an entry 
vehicle and to ascertain the effect of this precursor on the vehicle flow field. In order 
to accomplish this task, a computational method and program were developed to properly 
model the precursor. As part of this research, expressions were developed for the mass 
production rates of each species due to the photodissociation and photoionization 
reactions. Also, consideration was given to the absorption and emission of radiation and 
how it affects the energy in each of the energy modes of both the atomic and diatomic 
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species. Subsequently, a series of parametric studies were conducted covering a range 
of entry conditions in order to predict the effects of the precursor on the shock layer and 
the radiative heat transfer to the body. 

For this research, a nonequilibrium chemistry viscous shock layer (VSL) program 
developed for NASA 22 and modified extensively by Carlson and Gaily 23 was used to 
predict the conditions of the gas in the stagnation region behind the shock front. This 
program included nonequilibrium chemical reactions, multi-temperature effects, 
nonequilibrium nongray radiative transfer, as well as coupling between the radiation and 
the flow field. The radiation calculations also accounted for the fact that the excited 
electronic states of the atoms were not in local thermodynamic equilibrium. 

The radiative transfer calculations were carried out using a revision of a program 
called RADICAL, originally written by Nicolet 24 and later improved by individuals at 
the NASA Langley Research Center. However, it was extensively modified at TAMU 
to include nonequilibrium phenomena. Further, while the tangent slab approximation 
was used for the radiation calculations in RADICAL; as part of this research additional 
expressions were derived to correct for the geometric attenuation of the radiation in the 
precursor. RADICAL was also changed to include certain radiative processes which 
were present in the cold gas of the precursor but excluded in the original revision. The 
properties of the gas in the precursor were calculated using the one-dimensional Euler 
equations in conjunction with the radiation transfer calculations of RADICAL. For this 
study, the earth’s atmosphere was modeled as a nitrogen gas. These calculations 
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accounted for the production of each species by radiative reactions as well as thermal 
nonequilibrium between the electrons and heavy particles. 
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PRECURSOR FORMULATION 

As shown in Figure 1, the flow around a high velocity entry vehicle can be 
divided into several regions. Immediately behind the shock front, there is a highly 
nonequilibrium region dominated by collisional chemical reactions, in which the shocked 
gas relaxes toward to chemical and thermal equilibrium. Following the relaxation region 
are the near equilibrium and the equilibrium regions, where, as suggested by the names, 
the gas is close to or in thermal and chemical equilibrium. In these regions, the gas 
cools due to the emission of radiation. Finally, as the gas approaches the body, it enters 
the viscous or thermal boundary layer where, due to the cold vehicle surface, the gas 
temperature rapidly decreases with a subsequent increase in the density. 

A portion of the radiation emitted in the hot near equilibrium and equilibrium 
regions of the shock layer propagates into the region of cool gas ahead of the shock, 
where part of the radiation is absorbed. This radiation absorbed ahead of the shock 
wave, referred to as "trapped radiation", causes a heating as well as excitation, 
photoionization and photodissociation of the cool gas. This region of cool perturbed gas 
ahead of the shock wave is the precursor. Due to the precursor, the shock wave 
propagates into an excited, ionized and partially dissociated gas. 

Thus, the shock wave precursor is a nonequilibrium region dominated by the 
absorption of radiation and subsequent photoreactions. Flowever, with the relatively low 
temperatures in this region the collisional reactions have little influence on the gas in 
comparison to the photoreactions; likewise, the viscous effects in this region are 
negligible. As discussed previously, there was originally some confusion as to the 
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Figure 1: Nonequilibrium Shock Wave Structure 
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importance of diffusion in the precursor. Very near the shock wave, the effects of 
electron diffusion from the highly ionized gas behind the shock might be of importance; 
however, on the distance scale associated with the absorption of radiation and the 
photoionization and photodissociation processes, the effects of electron diffusion are 
negligible. 

The dominant absorption processes in the cool precursor ahead of a vehicle 
entering the earth’s atmosphere are the photoionization and photodissociation of 
molecular oxygen and nitrogen involving continuum radiation. Photoionization due to 
continuum radiation for atomic oxygen and nitrogen are also present to a lesser extent. 
The atomic photoprocesses are not as prevalent as the molecular continuum due to the 
relatively low population of atomic species compared to the molecular species. Likewise, 
atomic line absorption is also small compared to the continuum absorption. This behavior 
differs from several previous precursor studies which showed significant absorption 
through the atomic lines. 14,16,25 However, these studies were performed in monatomic 
gases in which the atom concentrations were high. 

In cool air, noticeable absorption of radiation begins at energies greater than 
approximately 6.7 eV. This absorption takes place in the Schumann-Runge band system 
of molecular oxygen 26 and the absorption cross section for this process is of the order 
10' 18 cm 2 27 . At 7.05 eV, continuum processes associated with the dissociation of 
oxygen become important and at frequencies above 8.55 eV, the Lyman-Birge-Hopfield 
band of molecular nitrogen absorbs radiation 28 . However, this molecular band 
corresponds to a forbidden transition, afTi g to X / E | + and is thus relatively weak, having 
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cross sections on the order of 10 ' 21 cm 2 29 . At frequencies above 9.76 eV, absorption 
occurs through continuum processes associated with the dissociation of molecular 
nitrogen through predissociation out of the excited dll g state 26 . This process will be 
discussed in further detail in the section on radiative processes. 

Photons with energies exceeding the ionization energy from the ground state of 
the oxygen and nitrogen molecules (/ 02 = 12.1 eV, I m = 15.59 eV) are also absorbed 
strongly in cool air. The absorption cross sections for these two processes are of the 
order 10' 18 cm 2 and change very little with frequency 27 . 

Gas Model 

For this preliminary study, the earth’s atmosphere was modeled as pure nitrogen 
rather than a nitrogen oxygen mixture. This approach is a common simplifying 
assumption when performing nonequilibrium, hypervelocity flow field calculations since 
the atmosphere is predominately nitrogen. In dealing with the precursor, however, the 
primary concern was whether or not the absorption processes of nitrogen sufficiently 
model those of air. After careful consideration it was decided that due to the 
predominance of nitrogen in the atmosphere it would be reasonable to represent the 
atmosphere as nitrogen in this initial study. 

The effects of thermal nonequilibrium in the precursor were included in this study 
by permitting the free electrons and heavy particles to have different temperatures. 
Further, it was assumed, that the free electrons and the electronic states were in 
equilibrium at a common temperature, which as discussed by Nelson and Goulard 15 , is 
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one of the limiting cases for the precursor. For this region of the gas, the temperature 
governing the electronic states would normally be expected to be greater than the heavy 
particle temperature but less than the electron temperature; thus, ideally a three 
temperature model should be used allowing a separate electronic temperature. 
Nevertheless, since the mechanisms and expressions for the transfer of energy between 
the electronic states and the free electrons are not known or well understood, it was 
decided to use only a two temperature model. However, in order to correct for the local 
thermodynamic nonequilibrium between the electrons and the electronic states, a collision 
limiting correction 30 was applied to the populations of the molecular electronic states 
when computing radiative emission and absorption phenomena. 

For this study, the mass production rates in the precursor zone due to collisional 
reactions were neglected in comparison to those due to photoreactions. The 
photoreactions used in this study included the dissociation of molecular nitrogen and the 
ionization of both molecular and atomic nitrogen, i.e. 

N 2 + h\ A 2 N 
N 2 + hv ** V2 N 2 + e~ 

N + hv N + + e~ 

The physical constants for the five species included in this study are given in Table 1 ; 
also, the elastic collisional terms in the electron/electronic energy equation were 
evaluated using the collisional cross sections of Gnoffo, Gupta and Shinn 31 . 

The effects of the absorption of radiation through free-free and bound-bound 
processes were also included in this study. While these processes do not cause chemical 
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reactions, they do cause an increase in the energy of the gas and their effects must be 
included in the electron/electronic energy equation. Absorption through atomic lines was 
neglected due to the expected low concentration of atomic species. 

Governing Equations 

As discussed above, the precursor is a region of the gas in both thermal and 
chemical nonequilibrium which is governed by the absorption of radiation, in which 
collisional chemistry and viscous effects are negligible compared to the radiative effects. 
The shock layer, on the other hand, is dominated by the radiative and viscous 
phenomena, as well as collisional chemistry effects. Due to this basic difference between 
the important aspects of these regions, it is efficient to calculate the flow field for each 
region separately. In this way, the viscous and collisional chemistry effects need not be 
calculated for the precursor where they are relatively unimportant. 

While the precursor is coupled strongly to the shock layer through radiation, the 
precursor is expected to affect the shock layer only slightly. In this study, the 
nonequilibrium radiation coupled shock layer flow field was first calculated assuming that 
the conditions ahead of the shock front were unperturbed. Then the resulting shock layer 
flow field was used to determine the radiative flux through the shock front and, using an 
iterative procedure, the precursor flow field as well as the radiation absorbed by the gas 
in the precursor. Finally, the precursor results were then used as freestream conditions 
for the calculation of a new shock layer solution to find the effects of the precursor on 
the shock layer. The methods used for the radiation and shock layer calculations will be 
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discussed in detail in other sections of the text. The remainder of this section will deal 
with the calculation techniques for the precursor region. 

The equations governing the fluid properties on the stagnation streamline in the 
precursor are the steady, one dimensional, nonequilibrium Euler equations. The basic 
one dimensional Euler equations are: 

Global Continuity 

4-(py) - o o) 

uX 


Momentum 

pV— + = 0 (2) 

etc dx 


Energy 

P V~ + = 0 (3) 

dx dx 

In equation (3), H is the total enthalpy of the gas defined in terms of the static enthalpy 
such that 

H = h + -V 2 (4) 
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The second term in equation (3) is the gradient of the radiative flux. This term accounts 
for the increase or decrease in the energy of the gas due to absorption and emission of 



14 


radiation. In addition to these equations, the equation of state for a two temperature gas 
is required, 


p = P j?r£pi— 
£ p 


D P - 

+ P-T7-— (T, - T) 
M e . p 


(6) 


To allow for the effects of thermal nonequilibrium, an electron/electronic energy equation 
was added to these equations, 
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In this equation, e,. is the kinetic energy of the free electrons, 3/2 kTJm e _, while e elcln and 
e° are the electronic and zero point energies of the nth species. The last three terms on 
the right hand side of equation (7) allow for the effects of the absorption of radiation. 
This equation is derived in detail in Appendix A. 
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Chemical nonequilibrium was accounted for in the precursor through the addition 
of a species continuity equation for each of the five species in the problem. These 
equations are of the form 

a 

pv- 

dx h\ dx 

where 



The term on the right hand side of equation (9) is the mass production rate of the nth 
species due to photoprocesses and is derived in Appendix B. The absorption coefficients, 
Ku K 2 and K 3 are those for the absorption and emission processes associated with 
each of the three photochemical reactions in the previous sections. 

In all of the above equations, the radiative terms, dq/dx, are the changes in the 
radiative flux due only to the absorption of radiation and not those due to the geometry 
of the problem. These terms will be discussed in further detail later in the text. 
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Solution Scheme 

The equations above are solved using a space marching technique starting at a 
point far from the shock front and marching in towards the shock wave. The point 
furthest from the shock wave was assumed to be far enough from the radiating shock 
layer that the ultraviolet radiation was absorbed between this point and the shock; thus, 
the gas properties at this point were forced to remain at the freestream conditions. The 
spacing between each spatial point was adjusted as the solution progressed so that the 
change in the flow properties between each point was greater than a preselected minimum 
value but less than a selected maximum value. This procedure forced a high 
concentration of points in the regions of large gradients and allowed the distance between 
points to increase in the regions of small gradients. 

At each individual point, the above equations were solved using an iterative 
process. For this iterative procedure, equations (1), (2) and (3) were integrated 
analytically from point, k-1, in the flow field to the next point, k, to give the following 
form: 

Continuity 
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Energy 

= <i3) 

As previously mentioned, the change in the radiative flux in the energy equation from 
point k-1 to point k is only that due to absorption and emission. 

The iterative scheme at each point, k, began by assuming that the flow properties 
are the same as those at the previous point, k-1. Then, the analytic expressions for the 
continuity, momentum and energy equations (11), (12) and (13) were used to calculate 
new values for the density, velocity and total enthalpy, respectively. Subsequently, the 
static/ total enthalpy relation, equation (4), was used to calculate a new value for the static 
enthalpy and the integral expressions for the species concentrations, equation (9), were 
integrated from the spatial point k-1 to the spatial point k using a fourth order Runge- 
Kutta (RK) integration scheme. Similarly, the electron/electronic energy equation was 
also integrated using the RK-4 scheme. At this point, the static enthalpy and 
electron/electronic energy expressions, equations (5) and (8), were solved for the heavy 
particle and electron/electronic temperatures, respectively. Finally, the equation of state 
. was solved for the new pressure. 

At this point, the change in each of the fluid properties between the new and the 
original values was calculated. If this change was greater than a selected convergence 
requirement, then the steps outlined above were repeated until the change in the 
properties through one iteration was within the convergence requirement. After the fluid 
properties at a point were calculated, the changes between the current point and the 
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previous point were calculated. If these changes were not within the maximum and 
minimum requirements selected then the spatial step size was increased or decreased, 
whichever was appropriate, by a factor of two and the process was repeated for this 
spatial point. If, on the other hand, the changes were within the maximum and minimum 
requirements, then the values for the current point were saved and the current spatial step 
size was used to move to the next point where the above process was repeated. 
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RADIATION FORMULATION 

In most of the previous work investigating shock wave precursors, there have 
been several assumptions imposed on the calculation of the radiative transfer. In some 
cases it was assumed that the shock layer emitted radiation as a black body at the 
equilibrium temperature behind the shock front . 14,25,32,33 However, the radiative flux 
passing through the shock front from the nonequilibrium shock layer is considerably less 
than that predicted by the black body assumption; therefore, the use of the black body 
assumption can lead to the prediction of larger changes in the precursor than are actually 
present. Similarly, several of the previous works used a multiple step absorption 
coefficient model 1,20,34 where at a given temperature the species radiative properties are 
assumed constant over specific frequency regions. Some researchers, such as Lasher and 
Wilson 3 , have ignored the spectral details entirely. However, since the precursor region 
is dominated by the photochemical reactions, variations in the radiative transfer might 
cause significant changes in the precursor. Likewise, the spectral details are very 
important in the precursor calculations since there are a number of important radiative 
processes which occur over different frequency ranges and the frequency of the photon 
absorbed as well as the process through which it is absorbed directly affects how the 
photons energy changes the energy of the gas. Without sufficient spectral detail, it is not 
possible to ascertain what portion of the radiation absorbed causes photoionization or 
photodissociation and what portion simply causes an increase in the internal energy of 
the gas. 
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Equations of Radiative Transfer 

Because of the necessity for accurate radiation predictions for the calculation of 
the precursor conditions, it was decided that a complete spectrally detailed method of 
calculating the radiative flux was in order. Thus, an extensively modified version of the 
program RADICAL was utilized. RADICAL, which was originally created as a method 
of predicting the radiative transport through the boundary layer of ablating entry 
vehicles, is also suitable for the calculation of the radiative flux through the entire shock 
layer. This program allows the user to select the frequency points to be used for the 
continuum radiation, so it was possible to obtain the spectral detail necessary for 
accuracy in the precursor calculations. RADICAL also performs detailed calculations 
of the atomic line radiation. 

RADICAL, like many of the schemes currently used in the calculation of radiative 
transfer, uses the tangent slab approximation. This assumption is a one-dimensional 
approximation of the full equation of radiative transfer. Basically, it treats the radiation 
emitted at a point in the gas as if it were emitted by an infinite plane of gas positioned 
perpendicular to the direction of travel of the radiation; this can be seen in Figure 2. By 
this assumption, the radiation emitted by the point, A, in the shock layer is treated as if 
it were emitted by the plane shown. As can be seen in this figure, this is a good 
approximation if we consider only the shock layer ahead of a blunt body. Since the 
thickness of the shock layer is much smaller than the body dimensions, each point in the 
shock layer is positioned close enough to the body that the rest of the gas in the radiating 
shock layer indeed appears to be of infinite extent. For the conditions of interest in this 
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study, the distance between the shock and the body on the stagnation streamline is on the 
order of 10 cm while the diameter of the body is 230 cm. 

The precursor, however, can extend to distances ahead of the shock on the same 
order of magnitude as the body diameter. For the conditions in this study, the precursor 
extends to distances on the order of 150 cm ahead of the shock. Therefore, in the 
precursor, the radiating shock layer no longer appears to be of infinite extent but instead 
appears to be a slab of finite diameter. 

In the one-dimensional problem, as in the shock layer, absorption is the only 
method by which the radiation is attenuated as it travels through the gas. Therefore, any 
decrease in the radiative intensity through the gas can be attributed to absorption, which 
in turn causes an increase in the energy of the gas equal to the decrease in the radiative 
energy. However, since the shock layer does not appear to be of infinite extent at each 
point in the precursor, the radiation no longer behaves one-dimensionally; consequently, 
in the precursor region the radiative transfer is a three-dimensional problem in which a 
decrease in the radiative intensity can occur because of the geometry of the problem as 
well as absorption. 

This effect can be better understood by considering a plane of radiating gas. In 
the one-dimensional problem this plane is of infinite extent and the radiative intensity will 
be constant at all points on a plane parallel to the radiating plane. However, in the case 
where the radiating plane is not of infinite extent, the radiative intensity on a plane 
parallel to the radiating disk will be a maximum on the centerline of the disk and will 
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decrease moving radially outward from this centerline. This attenuation of the radiation 
radially is due strictly to the geometry. 

Similarly, the radiation will attenuate faster in the direction perpendicular to a 
finite radiating disk than it would moving away from an infinite radiating slab. As shown 
in Figure 3a, the radiation emitted from an infinite slab travels only in the direction 
perpendicular to the slab. Therefore, since the radiative flux, q, is the quantity of 
radiative energy passing through a unit area, if there is no absorption the radiative flux 
will remain constant in the direction perpendicular to the radiating slab. This is a result 
of the fact that the radiative energy emitted in a unit area of the slab will continue to pass 
through the same area as it progresses perpendicular to the slab. However, for the case 
of a finite radiating slab there will be a decrease in the radiative flux moving in the 
direction perpendicular to the slab even though there is no absorption. This decrease is 
a result of the fact that the energy emitted from a finite slab will radiate evenly into the 
180 degree semi-sphere as shown in Figure 3b; therefore, as the energy emitted by this 
slab progresses outward the area through which it passes will increase, thus producing 
a decrease in the radiative flux. This decrease, however, since not due to absorption by 
the gas has no effect on the gas through which it passes. Therefore, in the case of a 
finite slab, the radiative flux will decrease due to both absorption and geometric 
attenuation. However, only the decrease due to absorption will cause an increase in the 
energy of the gas. 

In order to use RADICAL for the radiation calculations, it was necessary in the 
precursor to correct for the geometric attenuation of the radiation. The one-dimensional 
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Figure 3: Propagation of Radiative Energy from Finite and Infinite Slabs 
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equation of radiative transfer resulting from the tangent slab approximation is 


= 2 1 rC"E 3 (t < ) + 2nfs,E 2 (z v -z'jdz' - 2 

0 

( 14 ) 

T v 

- 2n fS v E 2 (x'-t,)*' 


In this equation, the first term is the radiative flux away from the wall, the second term 
is the radiative flux due to the gas between the wall and the current point of interest, the 
third term is the radiative flux from infinity and the fourth term is that due to the gas 
between the point of interest and infinity. 

If we assume that the point of interest is in the precursor and that the radiation 
emitted in the precursor is negligible, this equation reduces to 
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The first term here, as above, is that due to the radiative flux from the wall and the 
second term in this equation is due to the gas in the hot radiating shock layer. 

The equivalent equation assuming that the radiation emitted at each point is due 
to a finite diameter slab is 
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Comparing this equation and the one above, it can be seen that they are identical except 
for the terms in brackets which are dependent on the vehicle size and the location of the 
point of interest. These terms are essentially attenuation factors which correct each term 
in the tangent slab equation for the geometric attenuation associated with a finite slab. 
In this equation, the first bracketed term is the attenuation factor associated with the wall 
radiation and the second is that associated with the radiation from the shock layer. 

As discussed previously, the dominant absorption processes in the precursor are 
those causing photoionization and photodissociation of the gas, which occur in the 
ultraviolet frequency range. Since the radiation from the wall is in the infrared 
frequency range, it was expected to have little effect on the gas in the precursor. 
Therefore, for this study the shock layer attenuation factor was applied to all of the 
radiative terms. It is believed that this is a reasonable engineering approximation and 
should have minimal effect on the results of this study. Utilizing this attenuation factor, 
the radiative flux in the precursor can be expressed as 


<? v = ?v 


TS 


(17) 


where AF, is the attenuation factor and q, 75 is the radiative flux at the point using the 
tangent slab approximation. 

In the species continuity and energy equations, the terms involving the radiation 
appear as a divergence of the flux and are defined to account for the absorption and 
emission of radiation at a point. However, simple differentiation of equation (17) yields 

(IB) 


TS 


v , F 5 <7 v 

= AF.. — + q y 


dx 


dx 


ts dAF y 
' dx 



27 


In this expression, the first term on the right hand side is the change in the radiative flux 
due to the emission and absorption of radiation, but the second term is the change due 
to the geometry of the problem and should not affect the gas. Thus, it is neglected. 
Notice that if the second term was included in the species continuity and energy 
equations, an essentially transparent radiation would appear to be absorbed due to the 
spatial variation of the attenuation factor. 

A detailed derivation of the equations in this section is given in Appendix C. 

Radiative Processes 

In order to properly account for the effects of absorption and emission of radiation 
on the energy of a gas, it is necessary to have an understanding of how each radiative 
process physically changes the particles involved. The effects of the absorption and 
emission of radiative energy on the internal energy modes depends on the type of 
radiative process involved as well as the frequency of the photon absorbed or emitted. 
Radiative processes can be separated into three categories: free-free, bound-bound and 
bound-free. While the free-free and bound-bound processes cause an increase in the 
energy of the gas with no chemical change, the bound-free processes are associated with 
chemical reactions in the gas, such as photoionization or photodissociation. 

Free-free or Bremsstrahlung emission and absorption involves a free electron 
passing near a positive ion. This electron while passing through the coulomb field of the 
ion can emit radiative energy and through this processes suffer a decrease in kinetic 
energy. Similarly, the electron can absorb a photon of radiative energy which causes an 
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increase in its kinetic energy. It follows that the emission and absorption of radiative 
energy through the free- free process affects only the kinetic energy of the free electrons. 
Since there is no restriction on the kinetic energy which can be lost or gained by an 
electron, this is a continuous process and can produce photons at any frequency. 

Bound-bound emission and absorption is associated with a change in the internal 
energy state of the particle involved. In atoms, bound-bound emission and absorption 
is associated with a change in the electronic energy state of the atom; thus, emission and 
absorption of "line" radiation by atoms only affects the electronic energy of the gas. 
Since the electronic states are quantized and transitions between certain states are 
restricted due to quantum mechanical considerations, these bound-bound transitions only 
affect discrete frequencies. 

The emission and absorption of radiation due to bound-bound transitions in 
molecules is similar to that in atomic transitions. However, since molecules have three 
distinct modes of internal energy storage, electronic, vibrational and rotational, each or 
all can be affected by the emission and absorption of radiation. As with atoms, the 
absorption or emission of radiation by a molecule causes a change in the electronic state 
of the molecule. However, within each electronic state of a molecule there are a number 
of vibrational levels, as well as a number of rotational states within each vibrational 
level. Bound-bound molecular transitions involve a change from the original electronic, 
vibrational and rotational states to a new set of states and although the electronic state 
must always change, the vibrational and rotational states could feasibly remain 
unchanged. 
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Bound-bound molecular radiation is generally grouped into distinct bands. Each 
molecular band is associated with a transition between a pair of electronic states. 
However, these transitions can originate from any vibrational/rotational combination 
within the original electronic state and can end with certain restrictions governed by a 
set of quantum mechanical rules in almost any vibrational/rotational combination within 
the final electronic state. Since the internal energy states are quantized, the radiation 
associated with each pair of original and final vib rational/ rotational combinations appears 
at distinct frequencies and causes lines as in atomic bound-bound transitions. However, 
due to the large number of possible vibrational/rotational combinations within each 
electronic state, these lines tend to overlap and in most cases can be treated as a 
continuum process which occurs within a limited frequency range. Figure 4 is an energy 
diagram for the nitrogen molecule and shows the well known molecular bands of this 
molecule. 

The balance of energy in a molecule emitting or absorbing radiation through a 


molecular band is given by 


£ iow r, low _ law , 
elct + E vib + E rot 
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In this equation, the "low" values are the values for the lower energy state before the 
absorption of the photon and the "upp" values are those for the upper energy state after 
the absorption of the photon. If it is known through which molecular band the photon 
is absorbed, the difference between the final and initial electronic energies of the 
molecule, low , is known from knowledge of the two electronic states 

involved. The change in the individual vibrational and rotational energies can not easily 
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be determined without knowing the exact upper and lower energy levels of the two 
energy modes associated with the particular frequency of the absorbed photon. However, 
the sum of the changes in these two modes can easily be obtained by manipulating the 
above equation to give 

^ + - *v-A£. to (20) 

If the vibrational and rotational terms can be considered together, then the effects of 
absorption on these modes can be found using this relation. 

As previously mentioned, bound-free radiative processes are associated with 
photochemical reactions within the gas. As with free-free emission and absorption, the 
bound-free transitions are continuum processes. However, since the chemical reaction 
associated with a bound-free process requires a finite quantity of energy to occur, each 
bound-free process will have a lower limit on the frequency of photon which can be 
emitted or absorbed which corresponds to the energy required for the chemical reaction. 
Since the bound-free processes can emit or absorb photons of any frequency above this 
lower limit, they are continuum processes. Bound-free radiative processes can be divided 
into two groups depending on the type of chemical reaction involved, photoionization or 
photodissociation. The effects of photoionization on the gas will first be considered. 

While each species has one specific value for the ionization energy, the photon 
energy necessary for photoionization is dependent upon the electronic state of the particle 
absorbing the photon. Thus, the energy balance for a particle being ionized by the 

absorption of a photon is given by 

hv + E eU = I + ke e . 


( 21 ) 
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In this equation, hv is the energy of the photon absorbed, E elct is the electronic energy of 
the particle absorbing the photon, I is the ionization energy of this particle and ke e . is the 
kinetic energy of the electron freed in the ionization process. The limiting case which 
governs the minimum photon energy necessary for photoionization occurs when the freed 
electron has zero kinetic energy, ke e .= 0. If the particle in this case is in the ground 
electronic state, E elcl —Q, then the photon energy must be equal to the ionization energy 
of the particle. Further, if the particle is in an excited electronic state, then the photon 
energy required is less than the ionization energy by the energy of the electronic state of 
the particle and ionization from each successively higher electronic state requires 
decreasing photon energies. However, each successively higher electronic state is less 
populated; thus, the number of ionizations from the highly excited states is small. Also, 
as given by equation (21), if the photon energy is greater than that required the excess 
energy appears as kinetic energy in the departing electron. 

The photoionization process affects the distribution of energy through the various 
energy modes in several ways. Globally the energy of the gas increases by the photon 
energy, hv, however, the electronic energy of the gas decreases by the quantity Eelct due 
to the loss of the electronically excited particle. This particle splits into an ion and a free 
electron, which requires energy 7, accounted for in the zero-point energy of the created 
ion. Any excess energy in the process increases the kinetic energy of the free electrons 
through the kinetic energy of the created free electron. 

Photodissociation is a slightly more complex process than photoionization in that 
it involves the absorption of a photon by a molecule and the subsequent splitting of the 
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molecule into its atomic components. The discussion herein will deal only with the 
photodissociation of diatomic molecules since the gas involved in this study, nitrogen, 
forms diatomic molecules. Figure 5 shows the potential energy curves for the nitrogen 
molecule. 

The process of photodissociation is very similar to the absorption of radiation 
through the bound-bound molecular bands. As in the bound-bound process, the molecule 
absorbs a photon and transitions from its original electronic state to a higher state. 
However, if the photon absorbed is of sufficient energy the molecule will dissociate out 
of this upper state. If the upper electronic state in the dissociation transition has a 
potential well, as would be the case for a dissociation transition from the ground 
electronic state, X ! Y, g + , through the A + £ a + excited state, then the continuum radiation 
associated with this transition will adjoin the bound-bound molecular band associated with 
these two electronic states. On the other hand, if the upper electronic state in this 
process is continuous, or has no potential well, the continuum radiation will be 
independent of any molecular band. This situation would exist for the case of a 
transition from one of the lower electronic states of molecular nitrogen to either the 5 D* + 
or the 7 £ a + excited state. This is also the case if the minimum of the potential well on 
the upper electronic state is at a much greater intemuclear distance than that of the lower 
electronic state. However, the nitrogen molecule does not have a potential curve 
representing this case. 



34 



Figure 5: Potential Energy Curves for Molecular Nitrogen 
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The balance of energy between the initial molecule and the two atoms created in 


a photodissociation process is given by 
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In this equation, the "mole" values are the initial values for the molecule before the 
absorption of the photon. The energy of the photon absorbed is hv, the dissociation 
energy of the molecule is D, the electronic energies of the two atoms formed are 
E elc! am 1 and E etc “ m 2 and ke atoms is the combined kinetic energy of the two atoms. In this 
process, the total energy of the gas increases by the energy of the photon absorbed. 
However, the vibrational and rotational energy modes lose energy by the quantities E va , nu>lc 
and E ro[ molc , respectively; also, the effect of the photodissociation process on the electronic 
energy of the gas is given by 


» t- r atml j-atm2 ~ mole 

AE elct = E elct +E elct ~ E elct 


( 23 ) 


If it is known out of which electronic energy state the molecule dissociates, then this 
term can be evaluated since the electronic states of the atoms created are governed by the 
electronic state through which the molecule dissociates. The process of splitting the two 
atoms uses the quantity of energy D; this energy is accounted for in the zero point 
energies of the two atoms created. Finally, any excess energy above that necessary to 
cause the dissociation and excitation of the two atoms appears as kinetic energy of these 
atoms. 

Photodissociation of the nitrogen molecule generally occurs through a process 
called predissociation , 35 which is a radiationless process in which a molecule transitions 
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from a discreet electronic state to a dissociated state. The Franck-Condon principle 
dictates that for such a transition to occur, the potential curves for the discreet and 
dissociated states must intersect. Therefore, if a molecule is in a discreet electronic state 
which intersects a dissociated state and its vibrational level is great enough that the 
energy of the molecule is greater than the energy at the point where the two states 
intersect, then the molecule will transition into the dissociated state. 

Predissociation of the nitrogen molecule is known to exist from the dli g , B 3 Yl g 
and C 3 TL U excited states. The a 1 Yi g and B 3 H g states dissociate through transitions into the 
repulsive 5 '£ g + state while the C 3 U u state dissociates through a Tl u state which is not 
shown in Figure 5. 35 So, if a nitrogen molecule in the ground state absorbs a photon in 
the Lyman-Birge-Hopfield band with sufficient energy to place it in the dH g state above 
the point of intersection with the s E g + state, then the molecule will transition into the 
repulsive state and dissociate into two ground state nitrogen atoms. Likewise, similar 
phenomena occur for a molecule in the A 3 £ u + state absorbing via the 1st positive band 
to the B 3 n g state and transitioning to the repulsive 5 £ g + state or a molecule in the B 
3 II^ state absorbing through the 2nd positive band into the C 3 n u state and then 
dissociating through the state. All three of these dissociation processes produce two 
nitrogen atoms in the ground electronic state. 

For the relatively cool nitrogen present in the precursor, Herzberg 26 states that 
dissociation occurs primarily through the Lyman-Birge-Hopfield band and the subsequent 
predissociation out of the (^Ii g state. Therefore, for this study predissociation out of the 
and C J n„ states was neglected in an effort to reduce the numerical computations 
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required. It was noted during the course of this study that the absorption in the 1st and 
2nd positive bands which are involved in predissociation out of the B 3 Ti. g and C 3 II U states 
was negligible compared to that in the Lyman-Birge-Hopfield band. 

The radiative processes included in the calculation of the emission and absorption 
in the shock layer and precursor are given in Table 2. The radiative processes included 
in the shock layer are those which were originally accounted for in the modified version 
of RADICAL; these processes include not only the continuum processes, but also the 
atomic lines associated with the nitrogen atom. Since only continuum processes were 
included in the precursor, the continuum mechanisms originally included in RADICAL 
were retained; also, the photoionization of molecular nitrogen, the Lyman-Birge-Hopfield 
molecular band and the dissociation of molecular nitrogen through a continuum adjoining 
the Lyman-Birge-Hopfield band were added to the processes in RADICAL. 

The absorption coefficients for photoionization of molecular nitrogen and the 
Lyman-Birge-Hopfield molecular band were determined using theoretical expressions 
derived according to Zel’dovich and Raizer 36 . For the photoionization process, the 



where the photon energy, hv, is given in electron volts. The lower limit on the 
summation over the electronic states in this equation is governed by the requirement that 
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Table 2: Radiative Processes Included in the Shock Layer and Precursor 


Radiative Process 

Freouencv Ranee (eV> 

Shock Layer: 



Free-Free, Bremsstrahlung 

0.0 

^ he 

N - Low Frequency Ionization 

0.0 

£ he 

(Highly excited states) 

- High Frequency Ionization 

10.8 

<, he 

(Ground and first two excited 
states) 



- Atomic Lines 


— 

N 2 - Birge-Hopfield Molecular Band 

6.50 

r- 

c4 

VI 

2k 

X3 

VI 

- 1st Positive Molecular Band 

0.75 

< he < 4.5 

- 2nd Positive Molecular Band 

0.75 

< he < 4.5 

N 2 + - 1st Negative Molecular Band 

2.23 

< he < 4.46 

Precursor: 



Free-Free, Bremsstrahlung 

0.0 

<. he 

N - Low Frequency Ionization 

0.0 

^ he 

(Highly excited states) 

- High Frequency Ionization 

10.8 

<. he 

(Ground and first two excited 
states) 



N 2 - Ionization Continuum 

8.24 

<, he 

(Ground and first three excited states) 
- Birge-Hopfield Molecular Band 

6.5 

£ he < 12.77 

- 1st Positive Molecular Band 

0.75 

< he < 4.5 

- 2nd Positive Molecular Band 

0.75 

<, he < 4.5 

- Lyman-Birge-Hopfield Molecular Band 

4.77 

5 he sS 9.78 

- Dissociation Continuum 

9.78 

^ he 

(Adjoining Lyman-Birge-Hopfield 
molecular band) 



N 2 + - 1st Negative Molecular Band 

2.23 

£ he 4.46 
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the photon energy be greater than the binding energy for the state. Otherwise, the 
photon has insufficient energy to cause photoionization. 

For this study, the summation in equation (24) was limited to the lowest four 
electronic states of the nitrogen molecule; however, in the cool precursor the populations 
of all except the ground electronic state were small. It should be noted that equation (24) 
provides values near the ionization threshold on the same order of magnitude as those 
predicted by Zel’dovich and Raizer 27 as well as those predicted by Marr 29 . 

The absorption coefficient for the Lyman-Birge-Hopfield molecular band was 
found to be given by 

jy (113,314.97-n.6V0.U Jtv) 

k = 9.1458xZ0 16 e~ T ( 25 ) 

v T 

This equation was obtained from equation (5.113) of Zel’dovich and Raizer using an 
absorption oscillator strength of 3. 7x1 O' 6 from Allen 37 and then correcting to match 
experimental predictions given by Watanabe 28 . The absorption coefficient for the 
dissociation continuum adjoining this molecular band was assumed to be given by the 
expression 

= 4. 97x1 O' 20 N Ni (26) 

The constant in this equation was taken from the data presented by Watanabe for 
absorption through this process in cool air. 

Through the course of this research, several numerical studies were performed 
to determine the importance of various absorption processes in the precursor. For 
instance, originally predissociation was allowed through both the 1st and 2nd positive 
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molecular bands of nitrogen as discussed above. However, it was found through the 
course of this study that there was no significant absorption of radiation through these 
bands. Therefore, the effects of these dissociation processes were removed to decrease 
the computational effort required. 

Likewise, the Vegard-Kaplan band of molecular nitrogen was originally included 
in the precursor region for the frequency range from 2.45 eV to 5.84 eV. However, it 
was later found that there was insignificant absorption in this band; also, there was little 
evidence in the literature of the existence of this band in cool air. Similarly, although 
the Birge-Hopfield band was important in the high temperature shock layer, studies 
through the course of this work indicated that it was insignificant in the precursor; 
likewise, there was no reference in the literature to absorption through this band in cool 
air. It was also found that although photoionization was allowed for the ground as well 
as for the lower three electronic states of molecular nitrogen, ionization out of the ground 
state was much greater than that out of the excited states. Likewise, there was very little 
ionization out of the excited states of atomic nitrogen. 

Consequently, the dominant processes found in the precursor region of this gas 
were absorption through the Lyman-Birge-Hopfield molecular band, the adjoining 
dissociation continuum and the ground state ionization continuum of molecular nitrogen. 
There was no significant emission of radiation in the precursor and the molecular bands 
and ionization continuum associated with the excited electronic states were negligible. 
These observations are consistent with the predictions for the absorption spectrum of cool 
air from previous literature 26,27,28,29 
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SHOCK LAYER FORMULATION 

In order to properly model the precursor ahead of a shock wave, it is necessary 
to know the spectral details of the radiation which passes from the shock layer through 
the shock front and into the precursor. However, to calculate these spectral details the 
conditions of the gas in the shock layer must be known in detail. For the flight 
conditions of interest in this study, a number of important phenomena such as chemical 
and thermal nonequilibrium must be included to properly model the shock layer; also, 
since the effects of radiation are of primary importance in the precursor, it is desirable 
that they be included in the shock layer. The inclusion of these three phenomena can 
significantly affect the shock layer conditions which would in turn affect the radiation and 
hence the precursor. 

For this portion of the flow field, a viscous shock layer, VSL, scheme based on 
a highly modified version of the NASA code VSL3DNQ 22 was used. The VSL approach 
is basically a set of assumptions which when applied to the steady Navier-Stokes 
equations yield a more manageable set of equations which are valid in the region between 
the shock wave and the body. The benefit to using the viscous shock layer approach is 
that the VSL equations are solved by using a space marching technique starting at the 
stagnation streamline and moving aft around the body. This approach in general allows 
the equations to be solved faster than other techniques which require each cross flow 
plane, or the entire flow field, to be solved simultaneously. The VSL approach also 
alleviates the problems associated with solving the viscous and inviscid portions of the 
flow field separately and then matching the two solutions at the interface between them. 
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The viscous shock layer equations are commonly used to calculate the viscous flow 
conditions around blunt bodies and the case of interest in this study, an aerocapture 
vehicle entering the earth’s atmosphere from an outer planet such as Mars, involves a 
blunt body entering at hypersonic speeds. 

The VSL equations appropriate for the shock layer conditions in this study are as 
follows 22 : 

Continuity 



^-Momentum 



f 2 -Momentum 



( 29 ) 
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These equations are cast in a nonorthogonal, body-oriented coordinate system, Figure 6. 
The above equations were derived by nondimensionalizing the steady Navier-Stokes 
equations for a reacting gas and performing an order of magnitude analysis. In this 
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analysis, the normal velocity, u 2 , and coordinate, f 2 , were assumed to be of the order e 
where 


e 




\ p JJA 


N 


(33) 


and only terms in the Navier-Stokes equations of order e or less were retained while the 
higher order terms were neglected. 

The effects of thermal nonequilibrium in the shock layer were also included using 
a quasi-equilibrium electron energy equation 38 , 


1.23357s IQ' 10 
7*/sX[e-)m e 


vi> „ E, + w.—kT. 


e,eN I, 


N 


‘ 2 


(34) 


where 

SX = + (35) 

This equation was derived from the full free electron energy equation by Carlson and 
Gaily 38 by neglecting the derivative terms as well as the diffusion terms. The collisional 
reaction terms in this equation and in the species continuity equations were calculated 
using the collisional cross sections of Gnoffo, Gupta and Shinn 31 ; also, for the shock 
layer flow field calculations a lewis number, L e , of of one was used. 

The VSL program used in this study was modified extensively by Carlson and 
Gaily 23 . These modifications primarily involved the nonequilibrium chemistry and the 
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effects of thermal nonequilibrium. However, they also modified the code to allow the 
shock layer and radiation calculations to be coupled to the gas dynamics, thus 
incorporating the effects of the absorption and emission of radiation into the flow field 
solution. 
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GLOBAL SOLUTION SCHEME 

The computational methods used in the calculation of flow field conditions in the 
precursor, the flow field conditions in the shock layer and the radiative transfer were 
described in previous sections of this text. The means by which these three individual 
methods and programs were coupled together to calculate the entire flow field, shock 
layer and precursor, as well as the radiative transfer will now be discussed. The process 
of calculating the coupled flow field and radiative transfer was accomplished through a 
global iteration scheme involving all three of the programs discussed previously. This 
global iteration scheme is illustrated graphically in Figure 7. 

The solution process was initiated by first neglecting the effects of the absorption 
of radiation in the precursor. For these conditions, the VSL program was used to 
calculate the flow field conditions in the shock layer including nonequilibrium chemistry 
and radiative/gas dynamic coupling. This shock layer solution was then utilized to 
determine the radiative transfer in the precursor portion of the flow field. 

The radiative transfer in the shock layer and the precursor was first calculated 
with RADICAL using the shock layer solution and assuming that the flow field 
conditions in the precursor were constant at the freestream values. The absorption 
coefficients and radiative flux values from these radiation calculations were then used to 
compute the flow field conditions in the precursor region and this new precursor solution 
was compared to the previous flow field conditions. If the change between the two 
precursor solutions was greater than a preselected convergence criteria, the radiative 
transfer was calculated using the updated precursor flow field. This iteration scheme was 
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Figure 7: Block Diagram of the Global Iteration Scheme 
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repeated until the radiative transfer and the precursor flow field were fully coupled and 
converged. 

After calculating the coupled precursor flow field and radiative transfer, a new 
shock layer solution was calculated including the effects of the precursor ahead of the 
shock. These calculations were performed with the VSL program utilizing the precursor 
results directly in front of the shock as freestream conditions. This new shock layer 
solution was then compared to the original solution to determine the effects of the 
precursor on the shock layer. If this comparison showed significant changes in the shock 
layer flow field and the radiative transfer through the shock layer from the initial results, 
then the precursor flow field was calculated again using the new shock layer solution. 
This process, which was repeated until convergence, allowed a solution to be calculated 
for the entire flow field which fully coupled the precursor and the shock layer. 
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RESULTS 

This section outlines the results obtained for a total of six test cases at three 
altitudes, 72 Km, 75 Km and 80 Km. The atmospheric conditions for these three 
altitudes are given in Table 3. These cases are representative of "typical" conditions for 
an aerobrake vehicle entering the earth’s atmosphere upon returning from Mars. The 
freestream conditions for these six cases are given in Table 4. 

All of the results discussed in this section were obtained for the stagnation 
streamline of a 2.3 meter nose radius vehicle. The shock layer calculations were made 
using 52 points between the shock wave and the body and allowing for atomic local 
thermodynamic nonequilibrium as well as radiation/gasdynamic coupling. The radiation 
calculations were made using 74 continuum frequency points selected to provide good 
spectral detail in the ultraviolet absorption region of interest in the precursor. A wall 
temperature of 4^6 °K was used in both the shock layer and the radiation calculations. 
The precursor flow field was calculated using a convergence requirement of 0.05 percent 
and the changes allowed between each point in the precursor were confined to the range 
from 1 to 10 percent. 

The discussion of these results is divided into two sections. The first is a detailed 
discussion of case 1, 72 Km altitude and a freestream velocity of 16 Km/sec, which was 
the most extreme case for the precursor. The second section is a discussion of two 
parametric studies, the first over a range of altitudes and the second over a range of 
velocities. In that section, the variation of the precursor flow field properties with 
variations in the freestream velocity and altitude are discussed. 
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Table 3: Atmospheric Conditions for Three Altitudes; 72 Km, 75 Km and 80Km 


Altitude 

(Km) 

Temperature, T a 
(°K) 

Pressure, p w 
(dyn/cm 2 ) 

Density, 

(g/cm 3 ) 

72 


41.87 

6.659xl0’ 8 

75 

' 

25.75 

4.335xl0- 8 

80 


10.72 

1.999xl0‘ 8 


Table 4: Freestream Conditions, Shock Standoff Distances and Radiative Fluxes 


Case 

Velocity, V„ 
(Km/sec) 

Altitude 

(Km) 

Shock Standoff 
Distance, 

(cm) 

Radiative Flux 
Through Shock, 
(W/cm 2 ) 

1 

16 

72 

6.60 

1,385.0 

2 

16 

75 

6.72 

776.2 

3 

16 

80 

7.25 

264.5 

4 

14 

80 

8.69 

126.9 

5 

12 

80 

10.70 

65.9 

6 

10 

80 

11.14 

54.2 
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Detailed Discussion of Case 1 (Alt = 72 Km, V. = 16 Km/sec) 

Figures 8 to 14 show the heavy particle temperature, electron/electronic 
temperature, total enthalpy, electron/electronic energy, pressure, density and velocity 
variations through the precursor for case 1. Likewise, Figures 15 to 19 show the 
variations in the five species mass fractions through the precursor; also, the radiative flux 
through the shock front and the shock standoff distance for this case are given in Table 
4. The radiation emitted from the shock layer for case 1 was the greatest of all the cases 
considered. Thus, this case experienced the largest flow field perturbations in the 
precursor region. 

Figure 8 shows that the heavy particle temperature increased steadily through the 
precursor region as the shock front was approached from upstream. However, for this 
case the change in the heavy particle temperature was still less than one percent. On the 
other hand, there was a large increase in the electron/electronic temperature throughout 
the precursor region as shown in Figure 9. The large values for the electron/electronic 
temperature far from the shock were a result of photoionization in this region caused by 
photons at frequencies much greater than the ionization threshold; the large amount of 
excess energy in this photoionization process resulted in the creation of electrons with 
a very high kinetic energy. However, as can be seen in Figure 19, although the 
electron/electronic temperature was high far from the shock the electron mass fraction 
in this region was extremely small. 

The electron/electronic temperature increased steadily ahead of the shock to a 
maximum value of approximately 6,300 °K at a distance of 40 shock standoff distances 



Electron/Electronic Temp., T e /T M (-) 
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Figure 8: Heavy Particle Temperature in the Precursor Region; 
(Alt = 72 Km, V„ = 16 Km/sec) 



Figure 9: Electron/Electronic Temperature in the Precursor Region; 
(Alt = 72 Km, V„ = 16 Km/sec) 
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Figure 10: Total Enthalpy in the Precursor Region; 

(Alt = 72 Km, V a = 16 Km/sec) 



Figure 11: Electron/Electronic Energy in the Precursor Region; 
(Alt = 72 Km, V. = 1 6 Km/sec) 
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Figure 12: Pressure in the Precursor Region; 

(Alt = 72 Km, V a = 16 Km/sec) 
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Figure 13: Density in the Precursor Region; 

(Alt = 72 Km, V„ = 16 Km/sec) 
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Figure 14: Velocity in the Precursor Region; 

(Alt = 72 Km, V. = 16 Km/sec) 
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Figure 15: N 2 Mass Fraction in the Precursor Region; 
(Alt = 72 Km, V„ = 16 Km/sec) 
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Figure 16: N 2 + Mass Fraction in the Precursor Region; 

(Alt = 72 Km, V m — 16 Km/sec) 
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Figure 17: N Mass Fraction in the Precursor Region; 

(Alt = 72 Km, V a = 16 Km/sec) 
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Figure 18: N + Mass Fraction in the Precursor Region; 

(Alt = 72 Km, = 16 Km/sec) 
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Figure 19: e' Mass Fraction in the Precursor Region; 

(Alt = 72 Km, V a = 16 Km/sec) 
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ahead of the body. It then decreased rapidly to a value of 4,290 °K immediately ahead 
of the shock front. This decrease in the electron/electronic temperature was a result of 
the production of "low" energy electrons through photoionization caused by photons of 
frequencies only slightly larger than the ionization threshold. The production of these 
"low" energy electrons caused a decrease in the average energy per electron, hence a 
decrease in the electron/electronic temperature. That this decrease in the 
electron/electronic temperature was a result of the production of "low" energy electrons 
rather than due to a transfer of energy from the electrons through elastic collisions can 
be seen in Figure 11. If energy transfer through the elastic collisions were responsible, 
there would have been a decrease in the electron/electronic energy corresponding to the 
decrease in the electron/electronic temperature. The region of decrease in the 
electron/electronic temperature also coincided with the region of rapid increase in the 
number of electrons on the gas, Figure 19. 

The photons with energies near the ionization energy were absorbed rapidly in 
front of the shock since the strongest absorption region for a photoionization process is 
at frequencies near the ionization threshold. The higher energy photons in the weaker 
absorption range, far from the threshold, escaped to distances further from the shock 
where they were absorbed causing the creation of high energy electrons as discussed 
above. It should be noted that a similar decrease in the precursor electron temperature 
was also predicted by Foley and Clarke 19 for argon. 

Like the heavy particle temperature, the total enthalpy of the gas increased 
continuously towards the shock. The increase in the total enthalpy of the gas in the 



60 


precursor was less than 0.25 percent; thus, for this case the increase in the enthalpy of 
the gas due to the absorption of radiation was negligible. However, as can be seen in 
Figure 11, the change in the electron/electronic energy was not negligible. This energy 
increased from a value of essentially zero in the freestream to a value on the order of 
10 9 immediately ahead of the shock front. It should be noted at this point that 99 
percent of the radiative energy absorbed in the precursor affected the electron/electronic 
energy of the gas and only 1 percent of the energy affected the heavy particle 
translational, vibrational and rotational energies of the gas. Likewise, of the increase in 
the electron/electronic energy, 96 percent was involved with an increase in the zero point 
energy of the gas. Therefore, the majority of the energy absorbed in the precursor was 
involved with the ionization and dissociation of the gas. 

The pressure in the precursor, which was displayed in Figure 12, showed the 
same continuous increase towards the shock as did in the heavy particle temperature. 
However, due to the elevated kinetic energy of the electrons, the pressure increased by 
11 percent. Figures 13 and 14 show that in the precursor region, the density and 
velocity of the gas were essentially constant. This behavior verifies what was shown by 
Tiwari and Szema 20,34 and assumed by many other researchers. 13,14,1<5 - 33 

Considering the mass fractions of the five species it can be seen in Figures 15 
to 19 that the dominant chemical reaction in the precursor was the photoionization of 
molecular nitrogen. The mass fractions of the ionized nitrogen molecule was at least an 
order of magnitude greater than those for the nitrogen atom and the ionized nitrogen 
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atom; although, there were significant quantities of all three species immediately ahead 
of the shock. 

Due to the fact that the dominant change in the precursor region was caused by 
the photoionization of molecular nitrogen, the thickness of the precursor was considered 
to be the distance through which this reaction had an effect. By this definition, Figures 
16 and 19 show that for this case the shock precursor thickness was in the range of 75 
shock standoff distances, or 495 cm. Although there was a slight heating of the gas as 
well as the production of nitrogen atoms through photodissociation at greater distances 
from the shock, their effects were small compared to the changes within 495 cm of the 
shock. 

Figure 20 shows the frequency variation of the radiative flux at the shock front 
and in the precursor 495 cm ahead of the shock. The radiative flux at the latter point 
in the precursor is shown both corrected using the shock layer attenuation factor and 
uncorrected. Comparing the uncorrected flux at 495 cm in the precursor to the flux 
through the shock front shows the changes through this region due to the absorption of 
radiation; it can be seen that the dominant absorption in the precursor was at frequencies 
greater than the ionization threshold for molecular nitrogen, 15.59 eV. Although there 
was absorption at frequencies less than this threshold due to photodissociation of 
molecular nitrogen and photoionization of atomic nitrogen, the amount of energy 
absorbed in these processes was small compared to that absorbed in the photoionization 
of molecular nitrogen. This result agrees with the previous statement in discussion of 
the mass fractions that the dominant reaction was molecular ionization. 
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Comparing the radiative flux at 495 cm in the precursor that was corrected for 
geometric attenuation to that through the shock front shows that although only a small 
portion of the radiation was absorbed, very little radiation reached this point in the 
precursor. As can be seen from the figure, the geometric attenuation of the radiation 
dominated over absorption by the precursor gas. 

The radiation propagating through the shock wave from the shock layer into the 
precursor was distributed over a wide range of frequencies, Figure 20. As can be seen 
from this figure, a large portion of the radiative energy was in the infrared frequency 
range (hi> < 5 eV). Most of the radiation in this region was emitted by the entry body 
itself; although, embedded within the continuum radiation from the body were a number 
of atomic lines. Also, the peak of radiation near 3.5 eV was due to three molecular 
bands, the 1st negative band of N 2 + and the 1st and 2nd positive bands of N 2 . There was 
also a large quantity of radiative energy in the ultraviolet frequency range. The radiation 
in the frequency range from 8 to 10 eV was due to the atomic lines of nitrogen. That 
above 10 eV was due primarily to the Birge-Hopfield band of molecular nitrogen as well 
as the ionization continuum and lines of atomic nitrogen. Through the visible frequency 
ranges (5 eV < hv < 8 eV) there was very little radiative energy. 

Figures 21 to 29 show the heavy particle temperature, electron temperature, 
pressure, density and mass fractions in the shock layer for case 1 both neglecting and 
including the effects of the precursor region ahead of the shock front. As can be seen 
in these figures, the inclusion of the changes in the gas ahead of the shock layer had 
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Figure 21: Heavy Particle Temperature in the Shock Layer; 


(Alt = 72 Km, V. = 16 Km/sec) 
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Figure 22: Electron Temperature in the Shock Layer; 


(Alt = 72 Km, V, = 16 Km/sec) 
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Figure 23: Pressure in the Shock Layer; 

(Alt = 72 Km, = 16 Km/sec) 
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Figure 24: Density in the Shock Layer; 

(Alt = 72 Km, V, = 16 Km/sec) 





Figure 26: N 2 + Mass Fraction in the Shock Layer; 
(Alt = 72 Km, V m = 16 Km/sec) 
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Figure 28: N + Mass Fraction in the Shock Layer; 
(Alt = 72 Km, V*, = 16 Km/sec) 
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Figure 29: e' Mass Fraction in the Shock Layer; 
(Alt = 72 Km, V„ = 16 Km/sec) 
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negligible effect on the gas in the shock layer. Likewise, there was no measurable 
change in the radiative heat transfer to the body. 

The primary change due to the inclusion of the precursor was in the conditions 
of the gas immediately behind the shock wave. These effects can be seen more clearly 
in Figures 30 to 35 which show the electron temperature and the mass fractions in the 
chemically reacting nonequilibrium region. Neglecting the precursor, the mass fractions 
for the free electrons, ions and atoms were zero upon crossing the shock; however, 
including the effects of the precursor these mass fractions had nonzero values. Likewise, 
including the effects of the precursor resulted in a slight increase in the electron 
temperature in the region immediately behind the shock front. 

The case discussed above was representative of all of the cases calculated for this 
study and the remaining cases, being less extreme, had even less effect on the shock 
layer and radiative heat transfer to the body. Therefore, for the entry conditions of 
interest to an aerocapture vehicle entering the earth’s atmosphere upon return from Mars, 
the precursor region can reasonably be neglected in the calculation of the shock layer. 

Parametric Studies 

Figures 36 to 45 show the heavy particle temperature, electron/electronic 
temperature, total enthalpy, electron/electronic energy, pressure and mass fractions for 
cases 1, 2 and 3. All three of these cases were at a velocity of 16 Km/sec; however, 
each case was at a different altitude, 72 Km, 75 Km and 80 Km. The shock standoff 
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Figure 30: Electron Temperature in the Nonequilibrium Region of 
the Shock Layer; (Alt = 72 Km, V x = 16 Km/sec) 



Figure 31: N 2 Mass Fraction in the Nonequilibrium Region of 
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Figure 32: N 2 + Mass Fraction in the Nonequilibrium Region of 
the Shock Layer; (Alt = 72 Km, V a = 16 Km/sec) 
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Figure 33: N Mass Fraction in the Nonequilibrium Region of 
the Shock Layer; (Alt = 72 Km, V a = 16 Km/sec) 
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Figure 34: N + Mass Fraction in the Nonequilibrium Region of 
the Shock Layer; (Alt = 72 Km, V a = 16 Km/sec) 
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Figure 36: The Variation of Heavy Particle Temperature with Altitude in the 
Precursor Region; (Alt = 72, 75 and 80 Km, V a = 16 Km/sec) 


V„ = 16 Km/sec 



° ° 0 0 ° 72 Km 


— - 75 Km 

1 1 I 

<> o 80 Km 

i i t t i 


1 . 26 51 76 101 126 151 176 201 226 

Distance from Shock, X/X S hock (-) 

Figure 37: The Variation of Electron/Electronic Temperature with Altitude in 
the Precursor Region; (Alt = 72, 75 and 80 Km, V„ = 16 Km/sec) 
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Figure 38: The Variation of Total Enthalpy with Altitude in the Precursor 
Region; (Alt = 72, 75 and 80 Km, V m = 16 Km/sec) 



Figure 39: The Variation of Electron/Electronic Energy with Altitude in the 
Precursor Region; (Alt = 72, 75 and 80 Km, V a = 16 Km/sec) 
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Figure 40: The Variation of Pressure with Altitude in the Precursor Region; 
(Alt = 72, 75 and 80 Km, V*. = 16 Km/sec) 



Figure 41: The Variation of N 2 Mass Fraction with Altitude in the Precursor 
Region; (Alt = 72, 75 and 80 Km, V a = 16 Km/sec) 
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distance and radiative transfer through the shock wave for each of these cases are given 
in Table 4. 

From these figures, it can be seen that at a constant velocity as the altitude 
decreased the magnitude of the changes in the precursor increased. This inverse 
relationship between the magnitude of the changes in the precursor and altitude 
corresponded with the trends observed by Dobbins 25 as well as those observed by Tiwari 
and Szema 34 and was a result of two factors. First and foremost, as shown in Figure 46, 
with the decrease in altitude the radiative flux through the shock increased due to an 
increase in the extent of the equilibrium region in the shock layer. Second, with the 
increase in density at the lower altitudes, a larger percentage of the radiation passing 
through the shock was absorbed before being attenuated due to geometry. 

It should also be noted that as the altitude decreased, the length of the precursor 
region decreased. This change was a result of the increased density at lower altitudes, 
which caused the radiative mean free path to decrease. Hence, the radiation was 
absorbed in a shorter distance ahead of the shock. This trend was also predicted by 
previous studies. 15 ' 54 

Figures 47 to 56 give the heavy particle temperature, electron/electronic 
temperature, total enthalpy, electron/electronic energy, pressure and mass fractions for 
cases 3, 4, 5 and 6. All of these cases were at an altitude of 80 Km and the freestream 
velocities ranged from 10 to 16 Km/sec. From these figures, it can be seen that at a 
constant altitude, as the freestream velocity increased the magnitude of the changes in the 
precursor also increased. This trend was a result of the increase in the equilibrium 
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Figure 46: Variation of the Radiative Flux Through the Shock with Velocity 
and Altitude 
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Figure 47: The Variation of Heavy Particle Temperature with Velocity in the 
Precursor Region; (Alt = 80 Km, V. = 10, 12, 14 and 16 Km/sec) 



Figure 48: The Variation of Electron/Electronic Temperature with Velocity in 
the Precursor Region; (Alt = 80 Km, V„ = 10, 12, 14 and 
16 Km/sec) 
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Figure 49: The Variation of Total Enthalpy with Velocity in the Precursor 
Region; (Alt = 80 Km, V a = 10, 12, 14 and 16 Km/sec) 



Figure 50: The Variation of Electron/Electronic Energy with Velocity in the 

Precursor Region; (Alt = 80 Km, V a = 10, 12, 14 and 16 Km/sec) 
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Figure 51: The Variation of Pressure with Velocity in the Precursor Region; 
(Alt = 80 Km, V a = 10, 12, 14 and 16 Km/sec) 
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Figure 52: The Variation of N 2 Mass Fraction with Velocity in the Precursor 
Region; (Alt = 80 Km, V. = 10, 12, 14 and 16 Km/sec) 
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Figure 53: The Variation of N 2 + Mass Fraction with Velocity in the Precursor 
Region; (Alt = 80 Km, V^ = 10, 12, 14 and 16 Km/sec) 



Figure 54: The Variation N Mass Fraction with Velocity in the Precursor 
Region; (Alt = 80 Km, V a = 10, 12, 14 and 16 Km/sec) 
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Figure 55: The Variation of N + Mass Fraction with Velocity in the Precursor 
Region; (Alt = 80 Km, V a = 10, 12, 14 and 16 Km/sec) 



Figure 56: The Variation of e' Mass Fraction with Velocity in the Precursor 
Region; (Alt = 80 Km, V n = 10, 12, 14 and 16 Km/sec) 



85 


temperature in the shock layer as the velocity increased and the accompanying rise in the 
radiative flux through the shock front, shown in Figure 46; this trend is also in general 
agreement with the results and predictions of previous researchers. 13,25,34 It is also 
evident from these figures that the precursor thickness increased with velocity, again as 
a result of the increased radiative flux with velocity. As the radiative energy passing 
through the shock increased, a larger distance was required for this energy to be 
absorbed or attenuated ahead of the shock front. 

The results presented in this section give the same general trends as those shown 
in previous research for monatomic gases and for hydrogen. Likewise, the electron 
number densities immediately ahead of the shock front are of the same order of 
magnitude as those measured by Omura and Presley for a nitrogen gas. 12,13 In their 
research, they measured electron number densities of lxlO 12 /cm 3 for an 11.89 Km/sec 
case at a freestream pressure of 266 dyn/cm 2 in a nitrogen gas. The electron number 
density immediately ahead of the shock for case 5, p„ = 10.72 dyn/cm 2 and V a = 12 
Km/sec, was 5.5xl0 10 /cm 3 . The trend of the results presented herein with changing 
altitude and hence freestream pressure, show that the order of magnitude difference in 
the pressure for these two sets of results could account for the order of magnitude 
difference in the electron number densities. If the electron number density for case 5 is 
scaled linearly with pressure to match the pressure for Omura and Presley’s results, a 
value of 1.36xl0 12 is obtained. It should be noted, however, that Omura and Presley’s 
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results for air indicate that the presence of oxygen in the gas will have significant effects 
on the electron production in the precursor. Their results indicated that including the 
effects of absorption by oxygen in the precursor results in an order of magnitude increase 
in the electron number densities ahead of the shock wave. 
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CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE WORK 

A method has successfully been developed to calculate the chemical and thermal 
nonequilibrium precursor flow field as well as the subsequent effects of this region on 
the shock layer and radiative flux to an entry vehicle body. The precursor flow field has 
been calculated for a series of test cases representing flow field conditions for a vehicle 
entering the earth’s atmosphere upon returning from Mars. For these cases, it has been 
shown that the precursor region has negligible effect on the shock layer and radiative flux 
to the body. However, for an entry trajectory with a higher velocity or lower altitude 
than those considered in this study, the precursor could prove to be of significance. In 
any event, the effect of the hot radiating shock layer on the precursor is significant and 
leads to electron temperatures on the order of 6000°K and electron number densities as 
high as 7xl0 12 particles per cubic cm. 

It is suggested that research be continued in order to further improve the model 
used to calculate the precursor conditions in this study. Research by other authors has 
shown that the absorption of radiation by oxygen in the precursor can be significant and 
it is recommended that oxygen be added to this model in the future. The radiative 
processes of oxygen which are expected to have the greatest influence in the precursor 
are dissociation through the continuum associated with the Schumann-Runge band and 
molecular ionization. In addition, collisional chemistry involving election impact 
reactions may be important in the region immediately in front of the shock layer as also 
may be line radiation absorption. These phenomena should also be included in future 


studies. 
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APPENDIX A 

DERIVATION OF THE ELECTRON/ELECTRONIC ENERGY EQUATION 

WITH RADIATIVE EFFECTS 


For this derivation, the following notational scheme is used. The term, q, is the 
radiative flux vector. The negative of the change in this term, -dq/dx, is the total change 
in the energy of the gas per unit of time due to the absorption and emission of radiation. 
This change can be split into components, one for each individual energy mode of the 
gas, such that 


_3q = _fd£ N | 
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The terms on the right hand side of this equation represent the six possible modes of 
energy storage of the gas; free-electron kinetic, heavy particle translational, electronic, 
vibrational, rotational and the energy associated with the zero points of all species. 

However, at times it is also convenient to consider the change in the radiative flux 
due to each photoprocess. Therefore, the total change in the radiative flux is also 
considered by dividing into the changes due to each individual photoprocess, 
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The five terms on the right hand side of this equation are the changes in the radiative 
flux due to absorption or emission through the free-free processes, bound-free ionization, 
bound-free dissociation, bound-bound atomic lines and bound-bound molecular bands 
respectively. As shown, the sum of the changes due to each individual process is the 
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total change in the radiative flux. Equations (A-l) and (A-2) represent two methods of 
dividing the radiative energy absorbed or emitted by the gas. 

Utilizing the free electron energy equation derived by Carlson 39 and neglecting the 
diffusion, heat transfer and viscous terms since they are not being considered in this 
study, the electron energy equation including the effects of the absorption of radiation 
can be written as 
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The first term on the right in this equation is the increase in electron energy due to 
elastic collisions and the second term accounts for energy transfer due to inelastic 
collisions. The last two terms on the right hand side have been added to Carlson’s 
equation to account for the change in electron kinetic energy due to the absorption or 
emission of radiation and the transfer of energy between the free electrons and the 
electronic states, respectively. It should be noted that unlike Lee 40 the radiation term 
included above is not the total change in the radiative flux; as discussed in the 
"Radiative Processes" section of the text, the absorption and emission of radiation affects 
not only the free electron energy but all of the energy modes. Likewise, Lee does not 
include a term allowing energy transfer between the electrons and the electronic states. 
However, if it is assumed that the free electrons and electronic states are not in 
equilibrium with each other then a mechanism should be included by which energy can 


be transfered between them. 
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The equivalent equation for the electronic energy can be written as 


d_ 

dt 






(A-6) 


Now, for the precursor region it is assumed that the free electrons and the electronic 
energy modes are in equilibrium with each other and thus governed by the same 
temperature; therefore, it is appropriate to use a combination electron/electronic energy 
equation. This approach is also advantageous since it eliminates the energy transfer 
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terms between the electrons and the electronic states which are neither well understood 
nor well defined. Combining equations (A-5) and (A-6) leads to 
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Now, if a combined electron/electronic energy is defined as 



then this equation becomes 
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In the interest of computational efficiency, it is desirable to use the total change in the 
radiative flux rather than the individual components such as those in equation (A-9) since 
the evaluation of the two radiative terms in this equation will require a large number of 
frequency integrations. This effort can be decreased by using the total change in the 
radiative flux and subtracting off the portion which does not affect this equation. 

Since a portion of the radiative energy absorbed affects the energy stored in the 
zero points of the gas, adding these to the energy equation above is a step towards using 
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the entire change in the radiative flux. Multiplying the equation for the conservation of 
mass of each species by the zero point of the species and summing over all species yields 
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Adding this expression to equation (A-9) above and defining a new energy 
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leads to the following equation for the electron/electronic energy 
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By separating the mass production rates in this equation into those due to collisional 
reactions and those due to radiative reactions such that 

(A- 13) 

and then incorporating equation (A-l) for the radiative terms, equation (A-12) becomes 
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The remainder of this appendix involves simplifying equation (A-14) and deriving 
expressions for the last four terms on the right hand side of this equation. Before 
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proceeding with this derivation, a review is recommended of the section in the text, 
"Radiative Processes", which discusses the effects of each radiative process on the 
individual energy modes of the gas. A brief summery of the energy modes influenced 
by each radiative process is given by 


Radiative Process Energy Modes Affected 
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Defining 7/ = kf/k, lot for each individual radiative process rather than only for 
those causing chemical reactions, as in Appendix B, it is possible to find the change in 
the radiative flux due to each process 
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Likewise, since the energy of a photon is hv, the rate of photon absorption per unit 
volume per unit time for each process can be written as 
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These equations result from a derivation following the same steps demonstrated in 
Appendix B for the mass production rates. Equation (A- 16) will play a vital part in the 
development of the radiative terms in the electron/electronic energy equation. 

At this point, it is necessary to consider each radiative process and quantify its 
effects on each of the energy modes. First, the free-free process will be considered. As 
discussed previously, this process affects only the kinetic energy of the free electrons. 
By equation (A- 15), we have 
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Since the atomic lines involve only transitions between electronic states in an 
atom, the radiation emitted or absorbed through this process affects only the electronic 
energy. As with the free-free processes above, the effect can be determined with the use 
of equation (A- 15), 


= j dv (A- 18) 
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In this equation, the variable, Y BBA ' !ines , must include the absorption coefficients for all 
of the atomic lines. Equation (A- 15) was used to determine the effects on the energy of 
absorption or emission through the free-free and atomic line processes. Since each of 
these processes affects only one energy mode, equation (A-15) can be used to determine 
the amount of energy added to or taken away from the one energy mode involved. 
However, for the remaining processes it is necessary to find the effects of the absorption 
or emission of an individual photon on each energy mode and then use equation (A- 16) 
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to account for all of the photons involved. This is necessary since the energy absorbed 
or emitted through the remaining processes is distributed through several energy modes. 

•The absorption and emission of radiation through molecular bands involves not 
only the electronic states, but also the vibrational and rotational states. As mentioned in 
the text, the change in the electronic energy due to the absorption or emission of a 
photon through a molecular band is easily found from knowledge of the two electronic 
states associated with the band and is 

(A ' 19) 

where E upp elcti and E low elcti are the energies relative to the ground state of the molecule 
for the upper and lower electronic states associated with the ith molecular band. From 
this expression and equation (A- 16), the total change in the electronic energy due to 
molecular bands is 



The integration in this equation accoints for all photons absorbed or emitted in each 
molecular band and the summation is over all molecular bands in order to account for 
the effects of every band. 

The effect of absorption through the molecular bands on the vibrational and 
rotational energy of the gas is not easily calculated without exact knowledge of the upper 
and lower energy levels of the molecule associated with the absorption of the photon; 
however, these upper and lower states vary with the frequency of the absorbed photon 
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and are not easily found. The sum of the effect on these two energy modes can however 
be found easily. The change in these two energy modes is given by equation (20) as 
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As above the influence of all the molecular bands can be found with the use of equation 


(A-16) and by summing over the molecular bands, 
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The effects of the bound-free absorption processes on the zero point energy of the 
gas can be considered as a whole rather than splitting it into contributions due to 
ionization and dissociation. This is appropriate since the chemical reactions related to the 
bound-free photoprocess can all be represented by the reaction 
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The energy necessary for this reaction to occur is the difference between the zero point 
energies of the products and the reactants. For the ionization reaction this difference is 
equivalent to the ionization energy, 7, and for the dissociation reaction, it is the 
dissociation energy, D. For each photon absorbed through the ith process this reaction 
energy is given by 
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Using equation (A- 16), the increase in the energy associated with the zero points of the 
gas due to absorption by the ith bound-free process is given by 
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Summing over all of the bound-free processes results in the change in the zero point 
energy of the gas due to all of the radiative bound-free processes, i.e. 
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Now, substituting equation (B-17) for the mass production of each species due to 
photoreactions into equation (A-27), it is apparent that 
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The effects of the bound-free processes on the remaining energy modes will be 
found by first considering the ionization process and then the dissociation process. As 
discussed in the photoprocesses section of the text, the ionization process affects the 
electronic, free electron and zero point energies of the gas. For each photon absorbed, 
the electronic energy of the gas decreases by the electronic energy, E e icn of the particle 
absorbing the photon since the ion formed in the process is in the ground state and thus 
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has no electronic energy. Using equation (A- 16) and summing over the photoionization 
processes, the change in the electronic energy due to these reactions is 
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Due to the photoionization process, the energy of the free electrons increases by 
the kinetic energy of the electron created. This increase in the electron kinetic energy 
due to the absorption of a photon in the ith ionization continuum is given by equation 


(21) as 


- Av + E .i«, - h 


(A-30) 


Once again, the use of equation (A-16) along with a summation over the ionization 


processes leads to 


a ? f *> . > } a h 

rlr J. J h\3 dr 


(A-31) 


fa- i=1 0 


The final process to be considered is that of photodissociation. This radiative 
process affects the electronic, vibrational, rotational, translational and zero point energies 
of the gas. As discussed previously, the change in the electronic energy of the gas can 
be determined from knowledge of the molecular band through which the dissociation 
occurs; this change in the electronic energy due to a photon absorbed in the ith 
dissociation process is given by equation (23) as 


AC - zr atm 1 . r molc 

AC*/., + £, I , 


(A-32) 
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In this equation, the energy of the electronic state for each particle is considered relative 
to the ground state of the particle. Using this equation along with equation (A-16) and 


summing over the dissociation processes leads to the expression 

7 aan 1 ,-arm 2 r,molc\ 


*• up to n / Y?(lun l uOn J, _ 

m { h ' j 


dv 


(A-33) 


for the change in the electronic energy due to all of the bound-free dissociation 
processes. 

In a photodissociation process, the vibrational and rotational energy modes of the 
gas lose energy equal to the vibrational energy ,E molc va> , and cthe rotational energy, E 
mo!c ron of the molecule dissociated. It follows that the sum of the change in the vibrational 
and rotational energies of the gas due to the photodissociation process is given by 

X(-.C fc -0 a? fA . w 

UL + UL m l av & v 

As above, this equation results from utilizing equation (A-16) to account for all photons 
absorbed or emitted in each dissociation process and a summation to account for all 
dissociation processes. 

Lastly, the translational energy of the heavy particles is increased by the kinetic 
energy of the two atoms created in the dissociation process. From equations (22) and 
(23), the increase in the translational energy due to the absorption of a photon in the ith 
dissociation continuum is given by 


^ + E%‘ c + E™ k + hv — D : 


(A-35) 
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From this expression and equation (A- 16), the total change in the heavy particle 


translational energy due to the photodissociation processes is 

(dq)”** _ > } a q , 


dx 


Jtr 


■ £/ 
1=1 o 


h\ 


dx 


dv 


(A-36) 


The total change in the energy in each energy mode due to photoprocesses is 
given by adding the contributions from each type of photoprocess. This approach leads 
to the following expressions: 

Free Electron Kinetic Energy— 

1 dq^ 


, = /rf * * + £ 7 ^ ^ ^ (A-37) 

V^A,. J 0 dx - 


i= 1 


Av 




Heavy Particle Translational— 


(Id > J < ( - fa - g,) a g> (A . 38) 


Electronic- 


'i£'j 


i=l o 




S / 

1=1 n 


a* tr Av 

a? v 


ax 


/iv 0X 


Jv 


(A-39) 


fF ateI + r aIm2 _ r ,ru,/c 'j 

+ ^ / : — dV 

f=l J 0 


h\ 


dx 
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Vibrational and Rotational-- 


Wj + (M 

vib 




dx).,u \dx. 


rot 


Zero Point- 


ful 


hv dx 

-mole r-molc 


(A-40) 


n sF dia “ v p (- E,- E 1 P/7 
+ ^ J /v ' ^ £vt2>l d\ 

i=l 0 


hv 


dx 


d(j\ \—\ . ra^f o 

~ = -L w » 

K dx) e0 j=1 


(A-41) 


Substituting each of the appropriate terms into the electron energy equation (A-14) and 
simplifying leads to 

4(p«."> * + IX + <?. * *.-■£ 


af 


where 


etc 


n=l 


/» x 

E . coff o d — 
w n e„ - —-q 
n=l 5X 


(A-42) 
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(A-43) 
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APPENDIX B 

DERIVATION OF THE MASS PRODUCTION RATES 
DUE TO RADIATIVE REACTIONS 


As in Vincenti and Kruger 41 , the emission coefficient, j„ and the absorption 
coefficient, k„ are defined such that the rates of radiant energy emitted and absorbed per 
unit volume are given by: 


rate of energy emmited 
volume 


**> 4it 

/ fpj y dQdv 
0 0 


(B-la) 


rate of energy absorbed 
volume 


<» 4it 

f f pkJ^dQdv 
o o 


(B-lb) 


Since the energy of a photon is given by hv, it follows that the rate of emission and 


absorption of photons per unit volume is given by: 


rate of photons emitted 
volume 


“> 471 

f f—dQdv 
J J hv 

A A * * ' 


(B-2a) 


rate of photons absorbed 
volume 


If 


pVv 

hv 


dQdv 


(B-2b) 
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A set of radiative reactions can now be considered which are represented by the general 
radiative reaction 


n 


X 


T d ni x n + h\ * T a'X n (B-3) 

' ' n,i n X / njL n 

n=l / VJ n=l 

In this reaction, a n i and a / n i are the stoichiometric coefficients in the ith reaction for the 
species given by X n . Also, k ri and j ri are the absorption and emission coefficients for the 
photoprocess associated with the ith photoreaction. 

Using this notation, the total emission and absorption coefficients for the gas are 
given by 


n BF 




n 3F 


j\ = E>v i + h n 

M ‘ KBF ' 1 


(B-4) 


where y, nBF+; and k, nBF+1 represent the emission and absorption coefficients which are 
associated with the emission or absorption of radiation which does not cause a chemical 
reaction. This emission or absorption of radiation causes a change in the energy of the 
gas but no change in the composition. This is the case for radiation emitted or absorbed 
through the free-free process, the atomic lines or the molecular bands. Note that the 
summation in equation (B-4) is over the bound-free processes. 
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If the assumption is made that every photon absorbed by a given photoprocess 
causes the associated photoreaction to occur, then, using equations (B-2) and (B-3), the 
mass production of each species by the ith reaction is given by 


- ( a 'v- m .J/n£<*. 1 w. 1 >‘»* 


(B-5) 


where m n is the mass per particle of the species represented by X„ and is the density 
of the species absorbing the photon. The spectral integration in this equation is over all 
frequencies from zero to infinity. Looking at an individual reaction, the portion of the 
integral from zero to the threshold will have no effect on the mass production rate since 
the absorption and emission coefficients are equal to zero in this region. 

Now, under the quasi-equilibrium hypothesis the change in the radiative intensity, 


/„ for each radiative reaction is given by 


= -p^Vv-A,.) = 


(B-6) 


where 


k' r = k r (1 - e ~ h ’ ,kT ) 


(B-7) 


Using equations (B-5) and (B-6) the mass production rate of each species due to the ith 
reaction is given by 


a n) m nf f~r^~ 
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Now, if Y p ti is defined such that for each reaction, Y p wi = k' ri /k / rtot , then this equation 
can be expressed as 




SJdCldv 


(B-9) 


The superscript "p" is used on the ratio of the absorption coefficient for the individual 
process to the total absorption coefficient, Y r f , to indicate that this value is for the 
specific bound-free process associated with the photoreaction of interest. From equation 
(B-6) above, the total change in the radiative intensity due to all photoprocesses is given 
by 


= y 

& M 


dx 


n BF*\ 


= - E Kp,- s.) 

»'=l 

= -k' (L- SJ 


(B-10) 


Combining equations (B-9) and (B-10), the mass production of the nth species by the ith 
photoreaction is given by 


vv = -{a 1 - a ) m f j-^-dCldv 

n ‘ n J J /jv ^r 


(B-ll) 


0 0 


or, by rearranging the integral terms 



(B-12) 
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However, in this equation the integral over 4 t steradians of the change in the radiative 
intensity is the rate of change of the radiative energy flux, 


dq 4 r dl 

= f — -dQ 

J r)r 


(B-13) 


dx 


Applying equation (B-13) to equation (B-12) gives 


~ vP 


, Y : dq 

w = -(a .- a.)m f — - — -dv 
n ‘ n {hv dx 


(B-14) 


for the mass production rate due to the ith reaction. This equation gives the rate of mass 
production of the nth species due to the ith photoreaction. Thus, the total mass 
production of the nth species due to photoreactions is given by 


Tl BF 


w = w 

n / n. 
«=i 


■ - m »/ 


l BF 


E «,r “n,K 


o L i=1 


1 foy 
hv dx 


dv 


(B-15) 


or 


where 


■*. = ~ m n f 


KX, 

hv dx 


dv 


BF 


(B-16) 


(B-17) 


Y = ^(a' a )Y P 

i=l 

In this notational scheme, Y p r { is the fraction of the radiative energy emitted and 
absorbed at the frequency v due to the photoprocess associated with the ith photoreaction; 
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Y\ n , on the other hand, is the fraction of the radiative energy absorbed at the frequency 
v which is associated with the production of the nth species. The distinction between 
these terms is that Y p ri is defined for a specific radiative process and Y\„ is defined for 
a specific species. Defining Y\ n for each species rather than for each reaction is used 
only as a means of simplifying the equation since there may be many reactions associated 
with the production or destruction of each species. 

Equation (B-16) above is a general equation which allows the calculation of the 
mass production rate of each species due to photoreactions by knowing the change in the 
radiative flux due to absorption, the absorption coefficient for each individual process and 
the total absorption coefficient at each frequency v. We will now apply this to the 
specific set of photoreactions of interest in this study. From the precursor gas model 
described in the text, these reactions are: 

N 2 + h\ 2 N 

N 2 + h\ ** V2 N 2 + e " 

N + /zv ** V3 N* + e~ 



For this set of photoreactions, the mass production rate of each species is given by 
equation (B-16) above where 



(B-18) 
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APPENDIX C 


DERIVATION OF THE RADIATION ATTENUATION FACTORS 
In this appendix the equation for the radiative flux in the precursor region is first 
derived using the tangent slab approximation. This approximation treats the radiation 
emitted at each point as if it were emitted from an infinite plane of gas perpendicular to 
the direction of travel of the radiation. For the purpose of this derivation only, the 
radiation emitted in the cold precursor region is neglected compared to that of the hot 
shock layer. The equation for the radiative flux in the precursor is then be derived by 
treating the radiation emitted in the shock layer as if it were from a finite disk of gas. 
These two equations are then compared and attenuation factors obtained which modify 
the infinite slab equation for the geometric effects associated with a finite slab. 


The equation of radiative transfer in the s direction is given by 


— — 1r f 


= V('v-S v ) 


(C-l) 


where s is measured away from the point of interest and 7, is the radiative intensity 
towards the point. Placing both terms involving the radiative intensity on the left hand 
side and multiplying the integration factor 

(C-2) 
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Combining the terms on the left hand side of equation (C-3) results in the equation 

[ /.. e ~^ ds "\ = - kL S.. e ~^ ds " (C-4) 


d_ 

ds 


v“v 


In these equations, the s // terms are used to keep track of the various integration 
variables. Integrating both sides of equation (C-4) from the points to s and solving for 
the radiative intensity, I/s), results in 


l 3 ^ds" ('kids" s -(''kids" 

IJs) - /„(*,) e 1 - e 1 - J k[ s v e h ‘ ds’ (C-5) 


By combining the two exponential terms on the right hand side, this equation becomes 

(‘kids" s 

JSi v 


W -W e ' ’ - 




(C-6) 


To find the total radiative intensity at a point, the radiation due to the gas on both 
sides of the point is considered as shown in Figure 58. In this geometry the spatial 

variables s and x are related by the expression 

dx = -lids’ (C-7) 


where 

|i = Cos(d) (C-8) 

Applying this relation to equation (C-6) results in the following equation for the radiative 
intensity in terms of the spatial variable x 




w 




(C-9) 
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The radiative intensity at the point of interest can now be considered in two parts, that 
from the right, I r R , and that from the left, I, L . First consider the radiation from the right 
of the point of interest, between the point and the wall. With the geometry used in this 
derivation, the radiative intensity, 7„ is in the direction from right to left and the angle 

6 is in the range from 0 to 90 degrees. Also, Xj corresponds to the wall, 

x x = 0 


(C-10) 


W = 


wall 


Applying equation (C-9) to this portion of the radiative intensity leads to the expression 

-x % 

" /■ 


T \ T WCllt 

7 V (x) = 7 V e 


Jo n 


2 k's -f X i — dx" 

+ i-^e Jx » dx' 


(C-ll) 


for the radiation from the right of the point. The first term in this equation is the 
radiative intensity due to the wall and the second term is the radiative intensity due to gas 
between the point and the wall. 

Now, considering the portion of the radiative intensity from the left of the point, 
the x 1 and I,(Xj) values correspond to those at infinity, 

x i 

(C-12) 


w 


Applying equation (C-9) to this portion of the radiation yields 


-flid," x ic's 

/ v I W=/ > *e J -“ +[-Z±L e , -<‘ dx 

J I! 


(C-13) 
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for the radiative intensity from the left of the point. In this equation, the first term on 
the right is the radiative intensity from the boundary at infinity and the second term is 
that due to the gas between the point of interest and infinity. 

Introducing the optical depth, r,, defined such that 

= k[ dx (C-14) 


or 

t v = fk'dx' (C-15) 

0 

into the equations for the radiative intensity from the right and left results in the 
following expressions, 

_2» T » r frv-^v) 

7*(x) = l^e “ + f^-e » dx[ 

o ** 


<*y-Q T v s (*y~*y) 

/‘to = 7 v V “ + f-±e “ dz[ 

T u 


(C-16) 


(C-17) 


These equations for the radiative intensity are valid for both the infinite as well 
as the finite slab cases. The differences between these two assumptions occur in the 
integration of the radiative intensity over the solid angle to find the radiative flux. In the 
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tangent slab approximation this relationship between the radiative flux and the radiative 
intensity is of the form 

///.MM* (C-18) 

o -1 


or 


q™ = f fl v R \id\id( J> - j Jlf\id\id<b 

0 0 0 0 



(C-19) 


The angle </> in these equations is around the x-axis in the plane perpendicular to 
this axis. The radiative flux as defined in these equations is positive in the x-direction, 
or away from the wall. 

In equation (C-19), by integrating over the direction cosine, from 0 to -1 in 
the two terms, it is assumed that the radiation is emitted from an infinite slab, evident 
since a value of 0 for corresponds to a value of 90 degrees for the angle 6 shown in 
Figure 58. 

Utilizing equations (C-16), (C-17), and (C-19) and carrying out the integration 
with respect to </> from 0 to 2t, the radiative flux away from the wall for the tangent slab 
approximation can be written as 

l i T » (s— T v) 

q r = 2n f / v wa V* * d\i + 2nJ J S w e * dx' v d\L 

o oo 

-l _(vO -is (s-s) 

- In fl” \ie •* d\i - 2nJ J S v e * dx[d\i 


(C-20) 
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By switching the order of integration on the second and fourth terms and 
introducing the new variable 


H = -n 

d\i' = -d\i 

into the third and fourth terms, equation (C-20) can be written as 

1 1 JrX) 

q™ = 2rc/ v waU j H e * d\i +2 njS v fe * d\idx[ 


(C-21) 


i (<-0 


o o 


(C-22) 


- 2k/“ j* ^ + 2v: JS V je dp'dx'^ 


In this equation the integrals with respect to n are the 2nd and 3rd exponential integrals, 

l t_ 

E 2 (() = fe'^dv. 

° (C-23) 

1 _Jt 

£ 3 (f) = Jpe » d\ i 

o 

By introducing these expressions into equation (C-22) the equation for the radiative flux 
for the tangent slab approximation becomes 


?" = 2nl“° a E z (zJ + 2njs,E^-x[)dx[ 

0 

- 2K/ V "£ 3 (t"-T v ) + 2n J S v E 2 (x / v -x yi ) dx' y 


(C-24) 


This expression is the complete radiative flux equation for the tangent slab 
approximation. Up to this point no assumptions have been made regarding the point of 
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interest. However, by assuming that the point of interest is in the precursor, the 
integration in the second term of this equation can be divided into two parts, one over 
the shock layer and the second from the shock wave to the point in the precursor, 
resulting in the equation 

9,” = + 2 it J 

0 


* 2* / 5 v £ 2 (t v -x')Jt' - 2 7i 7“E 3 (t"-t v ) (C-25) 

+ 2 71 J S v £ 2 (t'-t v )Jt' 

m 

If it is also assumed that the radiation emitted in the cool precursor is negligible 
compared to that emitted in the shock layer, then equation (C-25) will further simplify 
to 



27cC J, £,(i: > ) + 2n 



(C-26) 


This equation gives the radiative flux in the precursor region due to emission in the shock 
layer by using the tangent slab approximation. However, as discussed previously in this 
text, this assumption is not valid in the precursor region. In the precursor, the shock 
layer no longer appears to be of infinite extent but appears to be a finite disk of gas. An 
equation will now be derived which properly accounts for the dimensions of the body and 
shock layer in the calculation of the radiative flux. 
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If it is assumed that the point of interest is in the precursor and that the emission 
from the precursor region is negligible, then equations (C-16) and (C-17) for the 
radiative intensity from the right and left side of the point reduce to 


c 

~ /* 


/* = I^e •* + 


/• 


(t v -c y ) 

dx' 


(C-27) 


I L = 0 


As can be seen in Figure 59, at a point in the precursor the shock comprises only 
a small part of the spherical field of view. Therefore, in this region the solid angle over 
which the radiative intensity should be integrated must be properly calculated. 
Considering the geometry of the problem, the expression relating the radiative flux to the 
radiative intensity is given by 

2 K 1 

= f f I*\id\id<b (C-28) 

0 

where /i bdy is the cosine of half of the angle subtended by the body, 

= Cos((5) (C-29) 

In equation (C-28) the radiative intensity from the left of the point of interest is omitted 
since this term is zero. By integrating over n in this equation from /i bdy to 1 rather than 
from 0 to 1 as in equation (C-18), the solid angle considered is only that subtended by 
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shock 



Figure 58: Geometry for the Derivation of the Radiative Flux from a 
Finite Diameter Body 
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the body. Combining equations (C-27) and (C-28) and performing the integration over 
<£, the equation for the radiative flux in the precursor region is 

4 v = 2 nl™ u f \ie » d\i + 2n j S v f e •* d\idx[ ( C " 3 °) 

0 Vbd, 


By dividing the integrals over jx into two parts, the first from 0 to 1 and the second from 
0 to as well as the introduction of the exponential integrals, this equation can be 
written as 


q v = 2 ic ir U E,(\) 


l*bdy 

J \ie 14 d\i 


1 - 




+ 2nfS v E 2 (x v -xi)dxi 


T v (t v * v ) 


1 - 


// 

0 0 


d\idx[ 


f E 2 (x v -x[)dx[ 


(C-31) 


L o J 

Now, introducing a change in the variable of integration for the two remaining integrals 
over n such that 


Vbdy 

d\i' = — — d\i 

V-bdy 


(C-32) 
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converts these integrals into the form of the exponential integrals 


9 V = 2izi: aU E 3 ^ v ) 1 - ^ 




W w,,' 




(*»-*») 


2nJs > £ 2 (T > -c , J )d< 1 






(C-33) 


« v = 2tt/”“£ 3 (x ) 1 
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r P (s- T v) ^ / 
t v Je 2 ^t v 
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/ E 2 (x^)dx'^ 


(C-34) 


Evaluating the integrals over r/ of the exponential integrals from 0 to r, s and 
substituting the definition for /z bdy , (C-29), equation (C-34) becomes 

. „v*U„ , J. „ 2 £ 3 K 5eC (P))l 


q v = 27c/ v w<u ‘£ 3 (t v )|1 - Co5 2 (P) 


* 3 CO 


(C-35) 
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This equation gives the radiative flux in the precursor due to emission in the shock layer 
properly accounting for the decreased field of view occupied by the hot shock layer. 
Comparing this equation to equation (C-26) for the radiative flux in the precursor using 
the tangent slab approximation, it can be seen that the presence of the terms in brackets 
is the only difference between them. These terms are essentially attenuation factors 
which can be used to modify each term in the tangent slab equation to account for the 
attenuation due to the three dimensional geometry of the problem. The first term in 
brackets is the attenuation factor for the wall radiation and the second term in brackets 




126 


VITA 

Scott Alan Stanley was bom in Shawnee Mission, Kansas on August 19, 1965. 
After graduating from Gainesville High School, Gainesville, Texas in 1983 he attended 
Cooke County College for one year. He then transferred to Texas A&M University 
where he received his BS degree in Aerospace Engineering in 1988. In the August of 
the same year, he returned to Texas A&M to continue his education in the Aerospace 
Engineering Department. He will receive his MS degree in December of 1990 and 
expects to begin a career in the Aerospace industry in the following January. 

Correspondence may be addressed to him care of W.C. and Nona Stanley, 709 
S. Morris, Gainesville, Texas 76240. 



